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i Abstract. We study asymptotic behaviour at time infinity of solutions close to the non- 

zero constant equilibrium for the Gross-Pitaevskii equation in two and three spatial di- 
mensions. We construct a class of global solutions with prescribed dispersive asymptotic 
behavior, which is given in terms of the linearized evolution. 
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1. Introduction 



Oh 

< 

• We consider the Gross-Pitaevskii equation for tp : M. L+d — > C 

id t ip = -AiP + (|0| 2 - (1.1) 

with the boundary condition 

\if)(t,x)\ — ► 1 as |x| — > oo. (1.2) 

This equation is a dynamical model for superfTuids and Bose-Einstein condensates, and 
has been extensively studied, especially concerning traveling wave solutions of the form 
if) = (f(x — ct), and dynamics of vortices (zeros of 0). See [2-10, 12-16, 20-25, 27, 28] 
and references therein. However we know very little about long-time dynamics of general 
solutions, for example about stability of vortices and traveling waves, or even of the constant 
solution ■0 = 1. Heuristically, the main difficulty is that small perturbations can continue 
^ | to interact with the non-zero background 1, and so do not easily disperse and decay. 

Thus we started in ^H] an investigation of large-time behavior of solutions tp = 1+ "small" 
as a first step toward understanding dispersive processes in this equation. The perturbation 
u = tp — 1 from the equilibrium satisfies the equation 

id t u + Au- 2Reu = F(u), F{u) := u 2 + 2\u\ 2 + \u\ 2 u. (1.3) 

The conserved energy and charge are written respectively 

.„ l2 (|u| 2 + 2Reu) 2 , f' 2 j 

|Vtt| 2 + — dx, / |u| 2 + 2Reu (fx. (1.4) 

For d = 2 and 3, unique global existence for the Cauchy problem has been proved first in 
P] for u G H 1 , and later in for any finite energy solution. 

In order to investigate dispersive properties of small solutions u, it is natural to linearize 
the equation around 0. The left hand side of ()1.3|) can be made complex linear by the 
following change of variable: 



V _:L «:=L7'~ 1 Re« + iIm«, U := y / -A(2 - A) -1 . (1.5) 
l 
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Then the new function v satisfies the equation 



id t v -Hv = -iV'hFiVv), H := V / -A(2 - A). (1.6) 

The linear evolution v = e~ iHt v(0) is expected to approximate small perturbations from 
the equilibrium. We have proved in the previous paper ^1 that this is indeed the case if 
d = 4 for small solutions v G H 1 . 

In this paper we turn to the physical dimensions d = 2, 3, where the dispersion becomes 
much weaker. Actually there exist traveling wave solutions with arbitrarily small energy 
for d = 2 [2], so it seems unlikely that the same result holds as in d = 4. However, it is still 
possible that small solutions u disperse in general if they are well localized in space at some 
time, since those traveling waves are spatially spread-out and do not belong to L 2 |T1|. Our 
theorems 11.11 11.21 show that at least there exist plenty of global dispersive solutions even 
for d = 2, 3. 

Since (|1.6|) is rather complicated, it seems natural to compare it with a simpler nonlinear 
Schrodinger equation (NLS) with the same nonlinearity: 

id t v + Av = F(v). (1.7) 

The scaling argument suggests that the quadratic nonlinear terms can be treated as a 
perturbation in L^°(L 2 ) only if d > 4. Hence we are led to work in weighted spaces, which 
provides more decay in time. Even with the optimal decay rate of the free evolution, the 
quadratic terms exhibit in general the critical decay order 1/t in l? x if d = 2. For d = 3, we 
can generally expect asymptotically free behavior for dispersive solutions. 

Coming back to our equation, the operator H has a singularity at frequency £ = 
similar to the wave equation, which is worse for time decay than the Schrodinger equa- 
tion. In addition, ()1.6|) apparently contains a singularity due to V in the nonlinearity. 
Nevertheless, it turns out that our equation is better than the NLS (|1.7|) . and in fact our 
argument for d = 2 does not apply to the latter, which appears simpler at first glance. This 
is because (|1.6|) has a special structure and thereby we can transform it to another equation 
with a derivative nonlinearity, in effect. We give the details of the transform in Section ri.il 

Before stating our results, we review the known results on the NLS for comparison. 
Dispersive global solutions have been constructed for the quadratic NLS in d = 2 only in 
the following two cases: 

((a)) F(u) = \ \u\u + Am 2 + X 2 u 2 , (\j G C), PUCE] 
((b)) F(u) = [Re(Au)] 2 , (A £ C), [IB] 

In the first case, the asymptotic profile is modified from the free evolution by a phase factor 
which depends only on Ao|u|u, while in the second case it is modified by the addition of 
a term with I? scaling concentration at £ = 0. Each argument depends essentially on the 
form of the modification, and so it seems difficult to combine these results to cover our 
F{u). For example when F(u) = |u| 2 , it is only known that there are no asymptotically 
free solutions in the L 2 sense with the natural decay property [2*§j . 

For d = 3, it is known ^7] that small initial data in certain weighted spaces lead to 



asymptotically free solutions in the above case (a) The final data problem (or construction 



of the wave operators) is easier and arbitrary quadratic terms can be treated in the same 
way as in Section |2 
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Now we state the main results in this paper. Hp and Hp denote the inhomogeneous 
and homogeneous Sobolev spaces respectively (cf. pQ), and we omit the subscript when 
p = 2. Denote by Bp q the homogeneous Besov spaces. In three dimensions, we have wave 
operators without size restriction. 

Theorem 1.1. Let d = 3 and < e be small (e < 3/68 is sufficient). For any T > 1, we 
define a Banach space Xj. by the following norm 

N*f = sups 1 / 2 - 8 !^^^!), (L8) 

where 1/p = We and 1/q = 1/3 — e. For any ip G H 1 satisfying \\e~ lHt (p\\x £ < oo, there 
exists a unique global solution i/j = 1 + u of Ql.l|) satisfying 

C(R; H 1 ) B V~ l u = e~ iHt ip + v', 

\W(t)\\(L rHi n W )(T,oo) < T- 1/4 - £ , \\v'\\ X e T < T-^l\ (L9) 

The above condition on ip is satisfied if cp G H 1 n H\, ^. We have the same result in the 
critical case e = if \\e~ lHt tp\\ x o is small enough. 

The threshold e = is related to the scaling property of the NLS with quadratic nonlin- 
ear ity in d = 3. 

In two dimensions, we can construct asymptotically free solutions for small final data. 

Theorem 1.2. Let d = 2 and if G H . Assume that 

(£}- 1/2 \^d k T^(0 G L°° n L? (1.10) 

for all multi-indices k > with \k\ = k\ + &2 < 2, and that ||y||^i is sufficiently small. 
Then there exists a unique global solution ip = 1 + u for (jl.lj) satisfying 

V- l u = z° -u + z' + z" G C(R; i/ 1 n if e ), 

^V, z' = t f e- iH ^\Uz°\ 2 ds, u = (2-A)- 1 U- 1 \u\ 2 , (1.11) 

•/ oo 

for any e G (0, 1) and t > 0, where the constants depend on e. 

Remark 1.3. The correction term ^ is coming from the normal form (see Sect. Il.lj) . It has 
a singularity at £ = 0, which can be worse than 1/|£|, because we do not know whether our 
solution u belongs to L 2 .. That is also the reason we describe it in terms of u, not ip. 

The correction term z' is essentially the same as in |18j for the NLS with (Reu) 2 , although 
we do not know whether it can be simplified as there, because of the singularity of our H(£) 
at £ = 0. It is probably not in L? x in general. 

However, these correction terms have no essential effect in the nonlinearity, and so they 
can be regarded as error terms if one does not require L? asymptotics for v or U2- 

In the next subsections, we explain our basic tools, namely the normal form and the L p 
decay estimate. 



z° = e- 
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1.1. Normal form. To eliminate the singularity at zero frequency £ = 0, we introduced in 
[16j the following transformation of normal form type: 

w = u+jP M!, (i.i2) 

where P was a Fourier multiplier cutting-off the higher frequency |£| > 1. The new function 
w satisfies the following equation 

iw = - Aw + 2Reu> + G(u), 

G{u) =(3 - P)u\ + Qui + PA\u\ 2 /2 + \u\ 2 u u (1.13) 

+ i[2Q(u x u 2 ) + VP • (n 2 Vni - uiV-u 2 ) + Q(\u\ 2 u 2 )}, 

where u = ui + iw 2 and Q = Id — P. It was crucial in |TB_ for d = 4 that Im G is essentially 
of derivative form. We also exploited the fact that the quadratic part does not contain u\ 
in the low frequency. 

Here we make a new observation that a special choice of P related to the equation 
leads to even better and much simpler nonlinearity. Let Q = U 2 = — A/(2 — A) and 
P = 1 - U 2 = 2/(2 - A). Then we have 2Q = -PA and 

G(u) = 2u\ + \u\ 2 ui - 2iV ■ P(u 1 Vu 2 ) + iQ{\u\ 2 u 2 ). (1.14) 

Hence the equation for z = V~ 1 w = U~ l {u\ + P\u\ 2 /2) + iu 2 is given by 

z = z° + f e- lH ^[N 2 (u) + N 3 {u)]ds, 



J oo 



(1.15) 



where we denote 

N 2 (u) := -2iu\ -2PU- 1 V- («iV« 2 ), iV 3 (n) := -i|u| 2 m + f/(|u| 2 n 2 ). (1.16) 

The new nonlinearity is roughly of the form (C/z) 2 + U(z 3 ). It is vital for our analysis in 
d = 2 that the quadratic terms consist only of derivatives. 

We will solve the above equation (|1.15|) for (z, u) and for t > T ^> 1 by the fixed point 
argument. Then solving the equation (|1.3|) for u from t = T by the result in ^U] (or by [2] 
for d = 3), and using local uniqueness of (#, «) satisfying 

z = e - iH{t - T) z(T) + f e- iH( t- s) [N 2 {u) + N 3 (u)]ds, u = Vz-P^, (1.17) 

we can deduce that our solution ^ := u + 1 satisfies the Gross-Pitaevskii equation (|1.1|) and 
extends globally in time. 

1.2. L p decay estimate. We recall the linear decay estimate proved in We call the 
pair of exponents (p, q) admissible if 2 < p, q < oo, (p,q) ^ (2, oo) and 2/p + d/q = d/2. 
We denote by q' = q/(q — 1) the Holder conjugate. 

Lemma 1.4. Let d>2. (i) Let 2 < q < oo and a = 1/2 — 1/q. Then we have 

\\e- UH <p\\ B o 2 < t- d °\\<p\\BO l 2 - (1.18) 
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(ii) Let (p,q), (pi,</i) and (p2,Q2) be admissible. Then we have 



-itH 



V\\ L vB0 o < C (p) IMIl2, 

q, 2 



-i(t-s)H 



f(s)ds 



U>1B° „ 

91 >2 



< C( Pl )C( P2 



(1.19) 



3^,2 



where C(p) is some positive continuous function of p, but diverges as p — > 2 w/ien d = 2. 

The above estimates are exactly the same as for the Schrodinger evolution e*' A . We had 
in ^B] some gain at £ = for d = 3, but we ignore it in this paper. The second last 
statement in Theorem 11.11 follows from the above estimate (i). 

For any s£i and T 6 1, we denote the full set of Strichartz norms of H s solutions for 
t > T by 



H\stz B T ■= sup C{p) \\u\\ L p( Tt00]H s). 
(p,q) -.admissible 



(1.20) 



When d = 3, this is just L°°H s nL 2 HI. When d = 2, it is slightly bigger than L^H'n^H^. 
We define the weighted Lebesgue space T by the following norm for any 0<6<l,sGM 
and T > 0: 



We denote the mixed norm by (where B is a Banach space) 

IMlLg, T (B) : = IIIK*)IIbJL. iT - 

The Holder inequality implies that 



J-Sl TS2 f- TS1+S2 

^6i,T A ^b 2 ,T L lj b 1 +b 2 ,T- 



We have also T C T , iff 



Moreover we have 



h < b 2 , s 1 + b 1 >s 2 + b 2 , T < T' . 



t- s £L s 0T (T>0). 



(1.21) 
(1.22) 
(1.23) 
(1.24) 
(1.25) 



The rest of this paper is organized as follows. In the Section 2, we deal with the three 
dimensional case, and the other sections are devoted to two dimensions. After explaining 
the main ideas in Section 3, we give the main bilinear estimate in Section 4, and then prove 
Theorem 11.21 in Section 5. 



2. Three dimensions 

In this section, we construct the wave operators in d = 3. The nonlinear terms are 
estimated simply by the Holder and Sobolev inequalities, and the wave operators are con- 
structed for the equation in our normal form by the standard fixed point theorem using 
the linear decay estimate. H 1 ! 2 regularity would be sufficient for the final state problem, 
but we do not pursue it in this paper. Once the solution u is constructed in C([T, do); H 1 ) 
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for some large T > 1, it is uniquely extended to a global one by the result in We will 
construct the asymptotically free solution by the fixed point theorem in the space 



(z,u) £ Stz^ n Xt, 



(2.1) 



for large T>1. 



2.1. The scaling critical case. We start with the simpler critical case e = 0. By using 
the L 3 decay estimate, we have for the quadratic term for t > T, 



1 \\e iH ^N 2 {u)\\ H ,ds < /* 1,-/1-' ^s-'lU' 



| s -t|-V^-i|| s V^( s )||^ lds 



<r II II 2 J.-1/2 



and for the cubic term 



^(s- i)Ar 3 (n)|| , g < 



S - t 



-1/2 -In 1/2 



u ( s )llH 1 H' U ( S )ll itx)rfs 



1 \\ u \\ yo 



<i- 3/4 IMIio 



X$U U \\L° 1/iT (L°-) £ J ll«ll^o||W||s tz i. 



The decay in S'tz is derived for the quadratic terms by 



JH(s-t) N 2^ di 



< \\N 2 



^3/4^(^3/2) 



< lit- 1 !!* 1 / 2 



and for the cubic terms by 



L°° 2,0 

3/4, T 



^ m— 1/2 11 ||2 n n 

< T I \\u\\ x0 \\u\\ su 



As for the normal form, we have 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



(2.6) 



for any p > 3/2. Thus we get the unique solution (z, u) for (|1.15j) by the standard fixed point 
argument in the space (|2.1jh provided that ||e~ t( p\\x° ls sufficiently small and T > 1. Then 
the solution is extended globally by the result in 0, and local uniqueness of the solution 
(z,u) in Ct{H l ) for p.!7jl . which follows easily from the Strichartz, Sobolev and Holder 
inequalities. 



2.2. Large data wave operators. Next we consider the case e > without size restric- 
tion. We define exponents q,q 2 ,q2 by 



l/q = 1/3 - e, 1/q 2 = 2/3 - 2e, l/q 2 = 1/3 + 2s. 



(2.7) 
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The decay estimate implies that for t > T, 



t 

e 



iH{t - s) N 2 (u)ds 



< / |t- s |-iAt*||JV 3 («(*))|| fll > 



H 1 Joo q2 

92 

l^ll - 1/2-8, 



(2.8) 



< T -i/2+2ei| 1 1 2 

rvj - 1 II"!! r l/i-oe/ H - 1 ^ • 



The Strichartz estimate implies that 

f e- m{t - s) N\u)ds < \\N\u)\\ L0 (H , 2) < r-^-IMI^ ( 2 . 9) 

The cubic term has additional T 1 / 4 decay due to the LfL 1 ^ bound by the same argument 
as in (|2.5|) . Then we use the complex interpolation 

[rl/2-2t/ ff h rl/4+EmUi _ r l/2-0/4-e(2-30) / ztIn /o ln \ 

where 6> € [0, 1] should be chosen to satisfy 

(l-0)/<7 2 + 0/6 = l/g, (2.11) 

i.e., 6* = 18e/(l + 12e) < 1. The last inequality is because e < 1/6. To embed the above 
space into Xf,, we need 

0/2 < We, 6>/4 - e(2 - 36>) > 2e, (2.12) 

i.e., 16e — 12e0 < 8 < 20e, which is satisfied by the above 9 with strict inequalities. In fact 
we have 

0/4-e(2-30) =2e + 5e/2 (2.13) 

Therefore we get T~ 5e ^ 2 as a small factor for the nonlinear term in the space (|2.1jl . The 
rest of proof is the same as in the critical case. □ 

3. Main ideas in Two dimensions 

In the rest of the paper, we deal with the case d = 2. In this section we describe the 
outline, and derive the key estimate in the next section, then finally prove the main theorem 
in the last section. 

3.1. Iteration scheme. Let u° := Vz° = Ve~ tHt (p. The integral equation is decomposed 
as follows 

z-z° = f e- iH ^N 3 (u)ds 

+ f e~ lH ^[N 2 (u) - N 2 (u°)]ds (3.1) 

J oo 

+ I e- lH ^- s) N 2 (u°)ds =: Tri(u) + Dif(u) + Asy(u°) 



The first two terms are estimated by simple Holder and Sobolev type inequalities, and the 
main task is to derive enough time decay for the last term Asy(u°), which is explicitly given 
by the data ip. In estimating Dif(u), we use i _1 L£° decay of u°, which forces us to assume 
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smallness of the data ip (this is usual in the case of critical decay). We further decompose 
Asy(u°) as follows: 

Asy{u°) = Asy'(u°) + z', z' = f ie~ iH{t ~ s) \Uz°\ 2 ds. (3.2) 

J oo 

z' is the only part where the oscillation of u° is completely canceled at £ = 0. 
3.2. Bilinear decay estimate. We will derive in the next section 

< ICr |fc| (0 1/2 Asy(u°) e r\\ogt) 2 H\ (3.3) 
Notice that we have by the simple t~ 1 L°° decay and the Holder that 

ip G H 1 H B®^ f e- lH{t ~ s) N 2 (u°)ds G L 00 ^ 1 ), (3.4) 

Jit 

so in l|3.3|) we are gaining roughly 1/t decay by losing £ at £ = 0, which is acceptable for 
our nonlinearity. 

The main idea of the decay estimate is as follows. For simplicity, consider the Schrodinger 
evolution H(^) = |£| 2 . Our quadratic terms are roughly of the form 

j* 9 \tm-v)${ri)dTid8 (3.5) 



in the Fourier space, where the phase function $ is given by one of 

$ := h(0 - H{t - v ) + H(n) = |£| 2 - le - v\ 2 + H\ 

*± := # (0 =F (£ - t?) + ^T(f7» = |£| 2 =F (|£ - r/| 2 + |r? 



$o corresponds to \u\ 2 and to u 2 and u 2 . We can gain 1/t by integration by parts in rj, 
picking up the divisor 1/|V^$|, where 

Vr,*o = -(ri-0 + V = Z, 

V,*± = =f[(»/-0 + »/]==F(2t/-0. ' 

Hence the singularity of 1/lV^ol i s canceled by |£| in (j3.5j) . We need to integrate twice, 
since we want to have t~ l after the integration in s. Then we get in the case of <3?o ; 

but it is almost in L 2 and so OK if we allow the loss of log t. 

In the case of <&±, depends on rj, which reflects the fact that u 2 and u 2 are oscillatory. 
But now we can integrate in s, because at the stationary point 77 = £/2, the phases <3?± do 
not vanish: 

= 2r/(£ - rj) = -2(r? - £/2) 2 - |£| 2 /2. (3.8) 

We are getting strong divisor 1/|£| 2 , but it is still OK around the stationary point, where 
we have |£| ~ |r/| ~ |£ — rj\ and the nonlinearity supplies \r]\\rj — £| decay. 

Since our actual symbol H(£) is degenerate at £ = 0, we get a stronger singularity, 
where the \r]\\ri — £| gain plays a crucial role. In addition, we should carefully compare the 
unbalanced radial and angular components. The detail starts in the next section. 
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4. Bilinear space-time phase estimates 



As seen above, we are going to have a non-stationary phase estimate for the bilinear 
expression with integration in space-time. Now we state the main estimate in a slightly- 
more general setting. 

Lemma 4.1. Let d = 2 and a > —1/2. Assume that F(^,ri) satisfies 



\d k F\ < 



\v\\£-v\ 



for all < | A; | < 2 and some nonnegative functions f,g£ L°° n L 2 and /',<7' £ L 2 , where 
m := min(|£ — r/|, |r/|). (/' and g' are not related to f and g.) Then we have 

r-t 



1^1 



e^ s F(^r))drids 



i 2 



(4.2) 



< r^iogtj^dl/lUconiall^lUcon^ + \\f'\\ L A\g'\\v) 



for < 9 < 1, t > 2, 



and 



^ = H(O + H( V )-H(n-O, 



<ii<e + {i + o-)±f- 



0-1 



< /i < 9+ (1 + min(fj,2o-))i±^ ($ = $±). 



(4.3) 



(4.4) 



The main part of proof is to derive precise lower bounds on the first derivative of the 
phase and compatible upper bounds for the higher derivatives. In doing that, we should 
carefully distinguish the radial and angular components, otherwise we would get too much 
singularity at £ = 0. 

4.1. Preliminaries. For any vectors £,r], we denote 

(0 = V± + W, [£} = V2TW, f=j|, Sr,--=Z-v, $ = Z-vSr l - (4.5) 

Then the phase function H(£) is written as 

H(t) = Hm = m}- (4.6) 

We will denote H'(£) := H'(\£\), etc. First we need to see that the above lemma applies to 
Asy(u°). Its Fourier transform is a linear combination of the form MG with 



G(S,r,) := 



'ivWv-m-'iv-zr 1 (forn 2 ), 

e • \V\(V - (for PU-'V • («iV«2)), 

^(£ - ■nWvi-n) (* = $o = (0 + - # (« - 0). 

J^(£ - r/)^) ($ = $+ = (£) - fffa) - H( V - £)), 
[^(-£ + r / )^(-r / ) ($ = $_ = #(£) + fa) + iffa - £)). 



(4.7) 
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Then in all six cases, our assumption (jl.lUj) implies the first condition of (|4,1|) with a = — 1/2 
and some nonnegative functions f,g€ L°° n L 2 determined by tp, and the second one follows 
from the assumption ip £ H . In addition, we observe that Asy'(u°) does not contain the 
terms with 3> = <J>q- By symmetry, we will mainly restrict our attention to the region where 



\v\ ^ \v ~ £l = m - 



(4. 



Since we are going to integrate by parts twice, we need up to the third derivatives of the 
phases. Let I(r) := H"{r)/r — H'(r)/r 2 . Explicit computations give us 



2(1 +r 2 ) 2r(3 + r 2 ) 12 
H (r) = — , H (r) = , H (r) = — ^, 



H""{r) 



60r 



r 



, J(r) 



4_ J/(r) _ 4(4 + 5r 2 ) 

"2 \ iy~\ 3 



(4.9) 



r 2 [r] J 



As for the differences, we have for any r > s > 0, 

H(r)-H(s) ~ (r)(r-s), H'(r)-H'{s) 



r 3 [r] 5 
r(r — s) 



ff"(r)-tf"(s)| < 
iJ'(r) if'fs 



< 



r — s 

TP 

r — s 
rs(s) 3 



|F'"(r)-F'"( S )| <— ^— ^(r-s), 
(r) (a) 

|/(r) -/(«)!< 



(r) 2 (a) 
r — s 
rs 2 (s) 3 



(4.10) 



For any vector v, we denote the partial derivative with respect to r] in the direction v by 
d v F(£,r]) := f • d v F(^,rj). We will omit the estimate with which is easier than that 
with <I> + . For the phases <l?o an d and for any vectors a, b, c, we have 

<S> = H(0±H( V )-H(r}-0, 
d a <S> = ±H'(i 1 )a ri -H , (r l -Oa v -t, 
d a d b $ = ±H"{r } )a T) b r) - H"( V - Os-cM 

± tm ± ± h'(t,-z) ± ± (4.ii) 

^ w »? I77 — £| 

d a d b d c $> = ±H'"(r])a ri b v c ri - H'"(r] - ^a^b^c^^ 

± I{v)( a , b, c) v - I(r) - £)(° 3 ^ c )»?-e' 

where the upper and lower signs correspond to $0 an d respectively, and (a, b, c) v denotes 
the symmetric 3-tensor defined by 



(4.12) 



We will use the following elementary geometry. For any r > s > and unit vectors a, /3, 
we have 



|ra — s/3\ ~ |r — s| + s|a — /3|, 



(4.13) 
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where < follows just by the triangle inequality and > follows by squaring the both sides. 
For any nonzero vectors a, b, we have 



|a + 6| 2 -(|a|-|6|) 2 (\a\ + |6|) 2 - \a - b\< 



\a + b\ 2 = 2(1 + a-b) 

mm mm 

Hence by putting (a, b) = (77, £ — 77) and (77, 77 — x), we have 

iei 2 -(N-|77-ei) 2 



(4.14) 



|r? - 77 -£| 2 = 2(1-77^) 
\v + rf^\ 2 = 2(1 + 7)^) 



(4.15) 



4.2. Estimate for |u| 2 , the case of <J>o- First we consider the phase $0, for which the 
integration in s does not play any role. For a fixed < 5 <C 1 and each £ 7^ 0, we split the 
integral region of 77 into the following three overlapping domains: 



A)(£) 



{77 € M d 
{77 G R d 
{77 G M d 



|77-e|>(l-2<5)|e|}, 
\r!-Z\\<(l-6)\{;\}, 
\t,-Z\<-(1-2S)\Z\} 



(4.16) 



and choose a partition of unity 1 = x+(v) + Xo( r ?) + X-(jl) satisfying suppx* C -D* for 
* = +, 0, — and 



KxM\<{H + \v-i\T k (W < 2). 



(4.17) 



Such functions can be given in the form x((M — \v ~ £l)/l£l) with some one-dimensional 
cut-off function Xi then its 77 derivatives are given by 



daX = 
d a dbX 



X J 



(0^ - a T? _g)(6 r? - 6^) „ 

| C |2 ^ + 



lei 



(4.18) 



The above bounds (|4.17|) follows from these identities, (|4.15|) . and |£| < I77I ~ (77 — £| in 

D±nD . 

By symmetry, it suffices to estimate only in D + and Dq. We first consider D + . Here we 
use the polar coordinates 77 = rO, or in other words, we choose the direction rj = 6 for the 
partial integration. By the definition of D + , we have 1 77 — £| < |?7| > |£|, and by (|4.15l) . 



Partial integration in r gives 
Fe im dr] 



(YF)e tm dr] = — \ (Y 2 F)e im dr] 



(4.19) 



(4.20) 
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where the operator Y is defined by 



i r 

Y:=-d r 



1 



+ 



r r d r $ rd r § (d r <S>) 2 <9 r <£»' 



Y 2 



8?<f> 



cl 3 $ | 2d r 



d^dr dl 



(d r $) 4 r(d r $) 3 (<9 r $) 3 r(<9 r $) 2 (d r $) 3 (d r $) 2 ' 
For the phase derivatives, we have the following estimates 
d r ^ = H'(n) - H'( V -0 + H'( V - 0(1 - 

- ${\v\ -\v-t\) + (v-${i- nn-6) Z . 



9 r 2 $ = F"(r ? )- J ff"(r ? -e) + 
|d 2 $ 



'0| < 



<»/> I 7 ? "CI 



\ri-Z\ 



|^-L |2 



<9 3 <l>o = tf» - ff'"(r, - + H'"( V - 0(1 - ^) 
+ 3/(r?-0^-?l^| 2 



l^ol < 



hi - I?? "CI , 1 - V« 



+ 



< 



(v-0 5 (v) Iv-tfiv-O 3 ~ 1^7 - el 2 ' 

where we have used |ry| > |ry — £|. Thus we obtain 

1 U+ >l ~ (6> T .^ )2| r? _ £|2 ^$0)2^ - + (^C,) 2 



< 



\v-£\(v) 



1-a 



\vW 2 \v-C\ 2 (0(ri-0 



(v) 



2-u 



\t\ 2 (t)\v\\v-t\(ri-t) 



T+^f{£-v)g(v)- 



(4.21) 



(4.22) 



(4.23) 



Next we consider D$. Here the main part of V<&o is hs angular component, which is not 
always close to either rj or rj — £, so we simply integrate in the direction of V^o- 1 Partial 
integration gives 



Fe tm dn = - / (AF)e tm dr] = — \ (A 2 F)e im drj, 



where the operator A is defined by 
AF := V- r^F 



V$ V$ • (2V 2 $ - A$I) • V$ 

VF i t^t- a '- F, 



|V$| 2 



|V$| 4 



(4.24) 



(4.25) 



On the other hand, our method for estimating in Do is not adequate in D+, where the angular difference 
terms such as \H'(r))/\r)\ - H'(r) - £)/\r) - £|| in gH) can not be controlled. 
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which satisfies 



11 ~ |V$| |V$| 2 ' 



|v$| 2 |v$p 



|V 2 $| 2 |V 3 $| 



(4.26) 



1^1- 



In the region D , we have \tj\ ~ |r/ — £| > |£|. Using (|4.13f) together with 1)4.91) and (|4.10j) . 

we have 

|v$ | ~ \h'(v) - H'( v - 01 + H'in - 01^-^1 

> (??)|^-r/-C|, 



|V 2 <J>o|<| J ff"(r / )-i/"(r ? -e)| + 



#'(77) i?'(r?-0 



+ 



#"(77) + 



\v\ 



\v\ I?? -CI 

I77 eg) 77 — ?? — ^ (8> 77 — £| 



|?7| 2 |t/| 

|V 3 d> | < |fT» - - 01 + |%) - I(V ~ 01 

+ (^"(^ + 1/(^)1)^-^1 



< 



+ 



\V\ 3 (V) ' |r/| 2 (r/) 3 
By the definition of Dq and (|4.15|) . we have 



|7? 77 51 Hh-ei m 2 



Hence we deduce that 



|V 2 $ | 



\v\ 



1 01 ~ |r/| 2 



and therefore 



\A 2 (xoF)\ < 



\v\\£-v\ x , \ 

|V^ | 2 |r/| 2 (0(e-r/) 1+,J (r ? ) <j/ ^ 



< 



iei 2 (0(^) 3+2, 

which is slightly better than the bound in (|4.23[) . 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



1-1 
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In conclusion, we obtain 
Fe i0ot dr] 



ri\>\t~v\ 



< 



f(£ - v)9(v) 



r/j/ 



+ 



+ 



i <i.h.-ci iei 2 (0(^) 1+2ff 77 

/(£ - »7)5(»7) 



< 



max(- log |£|,0) 



+ (0 



lei 2 

2-cr 



L°° \\9\\l°° + 



2/t\2+2<7 



l£l 2 (0 



where we used the Schwarz inequality and the condition a > —1/2 for the second integral. 
For 1 S> |£|, we have also 



/ 



\F\drj<f*g'. 



(4.32) 



Applying this estimate in the region {|£| < 1/t} and (|4.31f) in the rest, using the Young 
inequality L 1 * L 2 C L 2 , L 2 * L 2 C L°°, and appending the same estimate in the opposite 
region \rj — £| > |r/|, we obtain 



(0 1+CT kl / Fe Wot d V 



I 2 



< i- 2 (logt) 2 ||/|| LoonL2 || 5 || L ^ nL2 +i- 2 ||/'| 



L 2 \\9 \\L 2 , 



for t > 2. After integration in t, this estimate corresponds to the case 9 
remaining case < 9 < 1 is covered by interpolation, see Section 14.41 



(4.33) 
1 in The 



4.3. Estimate for u 2 , the case of <J> + . Next we consider the phase 3>+, for which we 
need to take account of the time oscillation, and so split the integral into more regions. Let 
X £ C°°(R) satisfy x( s ) = 1 for s < 1 and x( s ) = for s > 2, and denote x = 1 — X- Hence 
we have suppx C (— oo, 2] and suppx C [1, oo). We also denote 

C = r?-^/2, r = H, A=|r ? | + |r ? -C|-|el, 

(4.34) 

M = max(|7/|, |r/ — £|), m = min(|r/|, |r/ — £|). 

For a fixed positive i< 1, and each £ 7^ and t > 1, we introduce partitions of unity for 
77 G M rf by the following identities: 

1 = Xf + Xc, Xc = Xt "I" Xt + Xt + XX) Xx = xf + Xx +Xx + Xx, ( 4 -35) 
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and 

xf = x[ s-^j , xt = xc-x [s^x 



Xt = Xt- X±, Xt = Xt-Xo-X 



(0 2 A (4.36) 



Xx = Xx- x{t\C\), X X = XX- x{t\C\) ■ x(±4&, 

where x± and xo are the same as in (j4.16j) . Hence denoting by D* = supp xt > we have (see 
Figure HJ) 

l rf = D f UO+UD r u4uDf UD^UD X U L>^. (4.37) 
Remark that Dx and cover the stationary-phase regions in space and time, respectively. 




Figure 1. Decomposition of rj space for <J> + when |£| < 1 and t 3> 1. 

-D-r is an annular region, separating Dx into two connected components (if |£| < 1). The 
derivatives of the cut-off functions satisfy 

|VjxF|<M-l fe l l^-V^I + K^-V^l^m-*, 

iv.xxi < icr 1 + ^<€)ier 3/2 ^~ 1/2 , 

for | A; | < 2, where Z?t is identified with its characteristic function. We can easily derive 
these bounds using 

I 1 2 

m \v-Z\ m 2 (43Q) 

|V,A| = ^ + rPei~^, |VXTl<p£r72, 
together with similar estimates for the radial derivative in r\ — £. 
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In Dp, D^, and Dj,, we use the polar coordinates rj = r6 and integrate twice in the radial 
direction. 2 Since \rj\ > \rj — £| + |£|, we have 

-d r $ + = H'(r,) - h'( v - + H'(ri - £)(i + 



|<9 2 $ + | < + H»(r, - £) + ^ # l^-d 2 



< 



M + ^ ~ , 

fa) 1*7 — CI 



1 + 



|<£$ + | < H"'(r,) + H'"( V - £) + ~ Olfel 



< 



+ 



1 + Tfr-t 



In Dp, we have 



and also 



fa-O 5 \v-Z\ 2 (v-0 3 ' 

\rj — £\ 3> |£|, which implies that 

Hiei 



(0 



+ (v) > (v) 



rn 



(A; = 1,2). 



Therefore, defining Y by 1)4.20(1 with <I> = $ + , we get 



lY 2 ( F)l < £FMii7-ei/(g-*7)gfo) 

UF ~ (5 r $ + ) 2 |r ? -C| 2 (0(^-0 1+<J (^) CT 



/(£ - v)a('n) 



(0 



4+2<t 



Since a > —1, we obtain by the Schwarz inequality 



(0 2+max(<T,2,T) J XFFe^dq 



/ 2 



In D™, we have \rj\ ~ |£| > \rj — £| and by the radial component 



iei s 



to) (0 ' 

whereas in we have |"| ~ |£| ~ |™ — £| and by the angular component 

lei 2 



(0 



and we have (|4.42j) in both cases. Hence we have for * = 0, +, 



iei 4 i--ei 2 (0("-o 1+CT (r?r ~ iei 3 ir?-ei(^-o 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



(4.44) 



(4.45) 



(4.46) 



(4.47) 



2 Here we can not use the operator A because there is no cancellation between the angular terms H' ' (rj)/\r)\ 
and H'(rj — £,)/\r) — £j in V 2 <I?+, and also because the cut-off function has larger derivative in the angular 
direction. 
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For |£| > 1, this gives the same bound as in the last term in (|4,31|) . For |£| < 1, we use the 
shape of and the polar coordinates for r\ — £, 



\Y 2 (x* T F)\dv < 



< 



min(m,|g| 3 )||/|| L oo||ff|j 
|£| 3 m 
:°°||5'I|l°° log(l/|f|), 



-dm 



(4.48) 



where the factor min(m, |£| 3 ) is coming from integration in rj — £. We can treat D r in the 
same way by symmetry. Thus we obtain for * = 0, +, — , 



(0 



2 + (T / * 7-1 i$j.t J 



~ * 2 ||/llL 2 nL° c ll9llL2 nL oo. 



/ 2 



(4.49) 



Next we exploit the time oscillation in Z?x, which is split into the following two cases: 

fA»iei 3 /<o, 



a«|£| 3 /<0 2 , \\v\-\v-Z\\<(i-5M\- 



In the first case, we have \rf\ + \rj — £| ~ |£| < 1 and 



3 + )r? _ e |3 



lei 3 > iei 



fo] + V2 [ V -C] + V2 [£] + V2 ~ (0 ' 
while in the second case, we have |r/| ~ |r/ — £| ~ |£| and 

$+ = [H(0 - H(\ V \ + + \H{H + \V~ CI) - # w - ^ - 01 

ICI 3 



~ (0' 



where we used (|4.10p and the identity 

tti ,n m s tt / ~l\ o6(2a + 6) a6(a + 26) 
ff(a + 6) - (a) - (6) = - — v in , /, + 



[a + b} + [a] [a + b] + [6] 
Therefore we can integrate in s for £ 7^ 0: 



Af* := 



( X * X F)dsdr, 



for * = £?,+, 0, — . In particular, we have 

(0 ICI 2 



, , 1/ <i«i 3 (o 2+2 ° /( ^" )9( " )t! ' ? 



< 



\SK0 + h\ < i/t 

Hence by using the triangle inequality we obtain 

\\(0 2+2a x s \\ q < J 



ici < l/t 

<^ 3/2 H/llL-nL 2 ll5llL- 



1/2 



(0 

ici 1 / 2 



g\\L°°dQ 



(4.50) 



(4.51) 



(4.52) 



(0,6 >0). (4.53) 



(4.54) 



(4.55) 



(4.56) 
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In D^, we use the polar coordinate ( = le £ and partial integration in I: 



j Fe^dC = \ j (ZF)e^dC, 



where the operator Z is defined by 

The first derivative is given by 

d l <b + = H'(r ] )^ + H'{i 1 -0^l C 

= (H'( v ) - h'( v - 0)vc + H'(n - 0(^c + v^lc), 

where the radial component can be estimated by 

(H'( V ) - H'(rj - 0)VC > j|(M - \V ~ t\ fi + ff' C 

\ri\ 2 -\v-Z\H-t ^ \t\\<\% 



> 



(0 



\n\ ~ (0 ' 



and the angular component by 



2|Cl 



To rewrite it in terms of w, let a and /3 be the angles G [0, 7r] such that 

cos a = r] — £ • — £, cos (3 = rj ■ 
Then by the sine theorem and oo ~ sinw by \lo\ < ir/2, we have 



\v + V-£\ ^ sin a + sin /3 > 



1 



\v-£\ \v\ 



+ ^)|C||sinc| >^ 



m 



and so 



ff - h _ £)(%+ ,- 5<) >™^>«>^!. 



Thus we obtain 



> 



(0 : ™ 
The second derivative is estimated by 



m 



m 



w2 + Hei 



<0' 



~ (0 m 



lei , <™>M S 



(0 HCI 2 
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Then by using \r/\ ~ max(m, and the above estimate, we obtain 

< la^+Km^ + icr 1 ]. 

Thus we obtain 

l^(Xx F )l £ \*+\wx D xf(.Z ~ vMv), where 



(0 



III 



|£|m 



J x 



UJ 2 + 



4^ (0 14fl H l+ " 



11 Dt(Q 
m |C| rnVajf |3/2 



<0 1+<J l£l 2 <™) 1+<7 ICI (" 2 + $ 



1 + J_+ ^(0 



m |(| mV2|d 3 / 2 



Using |7/| ~ max(m, we have for |£| > 1, 
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(4.67) 



(4.68) 



+ 



(4.69) 



the first term is treated as the last term of Q4.31JI , and the second term by using the triangle 
and the Schwarz inequalities 



i>|<|>i/t 



i>ici>i/t Kh 



iS||/«/2-0IIl|IMIl- 



fd - vMv)d( 



Kl>l 



For |£| < 1, we have 



ICI 2 



/ i/(e-^)i 2 ^ 



, .ci>iKI 
^ ll/II^Uffll^. 



<4< 



+ 



■/?•? 



|C| 2 (^ + |C|2) | C |9/2| Cr 

and so in the polar coordinates £ = (i,u;), 



2n Idldu 



< 



i/tJo ^ 2 + iei 2 : 

log t + I log lell 



+ 



min(£, \^\' : 



-di 



Thus we obtain 



m 2 + ^n(aM) X+hm ^ m < t-l(logt)2||/|| LoonL2 || 5 || ioonL2 . 



(4.70) 



(4.71) 



(4.72) 



(4.73) 



We have the same bound for X by symmetry. 
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In Dj{, we use the operator A defined by Q4.24JI with $ = <J? + . Here we have |£| ~ 1771 r 
|^ — 77| and |Q| < 1/2, which implies by the same argument as in Q4.63[l that 



|r7+rPC|>|}. 



Hence we have 



\Vr,* + \>H'(ri-Q\Ti+Ti-Z\ > 



(OKI 

lei ' 



I „ +1 ~ |e| ~ |c| 



and so 

(0 lei ici 2 



w x . 



Id 3 (OKI (0 



2+2cr 



1 , (QDt 
ICI l£l 2 



< 



(0 2+2<T ICI 



1 , (Pgr 

ICI lei 2 



This bound is better than (|4.68[) . and so in the same way we obtain 

\\(0 2+2rT x°\\q m >i/t) < t-HiogmfU^yu^. 

On the other hand, we have for |£| < 1/t < 1, 

lei 21 



JJ \ X xF\dsd V 



< 
<t- 3 
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v\ < ICI 

|M|l° 



L°° \\g\\L°°drids 



(4.74) 



(4.75) 



(4.76) 



(4.77) 



(4.78) 



Putting the pieces together, we obtain 



(0 



2+min(<r,2o-) 



Fdrjds 



<t 1 (logt) 2 ||/|| L oo ni 2||5r|| i0 c ni 2. 



(4.79) 



Here the condition a > — 1/2 is inevitable for $ + . The difference from the case of $0 is the 
inbalance between the |£| from U and the 1/|C| from the partial integration, which costs 
one regularity. On the other hand, <I> + is better at £ = than <I>o, which is non-oscillatory 
at £ = 0. 

4.4. Interpolating estimates. Thus we have proved Q4.2|) in the case 9 = 1. On the other 
hand, we have for any k E Z, 



/ 



\F(C,rj)\d V 



< 2 k 



L 2 (2 k <\^\<2 k + 1 ) 

which implies that for any t > 0, 



\F{s,ri)\d<n 



<2 fc ||/'|| L2 || 5 '|| i2 . 



(4.80) 



<t\\f'\\ L 2\\g'\\ L 2. (4.81) 

2i ^ L2( 2 fe < | € | <2 fe+i) 

This corresponds to the case 9 = —1 in (|4.2|) . although it is divergent for the integral i — > 00. 
By applying real interpolation or the Holder inequality to dyadic sequences on |£| and t, we 
obtain the desired estimate Q4.2|) for < 9 < 1. □ 
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5. Proof of the 2D theorem 
5.1. Bootstrap setting. We decompose our solution z = z° + z 1 and u = u° + u \ where 

z°:=e- im <p, u°:=Vz°, (5.1) 
and starting from the above estimates, we will derive 

W^Wm Z *~ a , VW&w Z t~ , \Hu < t~ p , (5.2) 

with some a, (3 satisfying 

< 1/2, 1 - < a < 2/3, (5.3) 

by the standard iteration argument, (a, f3) can be arbitrarily close, but not equal, to (1, 1/2). 
We fix k £ (0, 1/4) such that 

l/2 + K<a<2/3- k. (5.4) 

Let H^ll^i < <S < 1. The L p decay Q1.18|) implies the following bounds on the free part 
u = u® + iu® '■ 

h (t)\\ m < y\\ H i, + ||v« ||l- < t-'yw^ < sr 1 , 

ll« (t)llL4<*- 1/2 ||^|lB0 A2 - (5 ' 5) 

The last quantity is finite for high frequency by interpolation of B\\ and H . The low 
frequency part is also finite, because ljl.lfl|) implies that 

(V>-^e ASSET 8 n^SST 1 , (5.6) 

for all p > 1 and q > 1 , which is proved simply by estimating the inverse Fourier transform. 
In the following three subsections, we derive estimates on the normal form, the trilinear 
terms and the quadratic difference terms, where we need not assume that u is the solution. 
For any function u, we denote 

\Hz> T HMU L 4 + II R*u\\ La Ll + ||v«|U 5 Ll , 

(5 7) 

\W\zt : =II^II l ^ L 4 + II B*u°\\ L i tTLs> + \\Vu°\\ l 1 itLs> + \\u - u°\\ z ^. 

Remark that is not a norm, but it is designed to measure different types of decay of u° 
and u — u°, namely dispersive and dissipative. Since (3 < 1/2, we have 

IMIl^ ~ W u Wzt- (5.8) 



5.2. Normal form. The quadratic part is estimated just by the Holder inequality: 

(5.9) 



L,|2|| <r* rrt— K|L,||2 <^ rp — K 1 1 „ , 1 1 2 

\U\ \\tol t2 >, 1 \\U\\ r a tJ 7 , 



0,T v 

|2 |„..|2| 



M ~H IU« r £3 £ T K (\\u\\z T + \\w\\z T )\\u-w\\ z > r . 
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5.3. Trilinear term. For 1 < p < oo, we have 

II^MIIhi £ IIMVII^ + \\u 2 Vu\\ LP . (5.10) 

We apply this estimate after expanding u = vP + u 1 , choosing different p for each term. For 
example, (ii2) 2 Vu2 is expanded into the following spaces 

||™Vu°|| L?/U? < ||u 2 ||% J|V^|| l i tL «,, 



1,T x 0,T sc 

«§«§Vt4|| £ a+l/3 L2 < ll^ll^ L JV«§||^ 3 TL? ||V^|| £ a T ^, (5.H) 



0,T x 

\\u 2 ulVul\\ 2a +0-3/4 r 4/3 < ||u 2 |L/3 L4 ||V«2||ia L 2 , 

3/4, T 31 0,T x 0,T x 

where we used the interpolation inequalities of Gagliardo-Nirenberg type: 

^ Ihlli^llVull^, \\u\\ L 8 < Ikll^llVnH^ 2 . (5.12) 
The other terms containing u\ are estimated in the same way. Thus we obtain 

|3 



\\Tn{u)\\ stz}r <T-^\\u\\% T , 
\\Tri{u)- Tri(w)\\ stz i <T- a -«[\\u\\ ZT + \\w\\ ZT } 2 \\u-w\\ z , 



(5.13) 



u = u° + u 1 



5.4. Quadratic error term. The quadratic difference term can be expanded by putting 

(5.14) 



= -2i{2u\ + u\)u\ - 2P[/~ 1 V ■ {u\Vu\ + u\Vu° 2 + u\Vu\), 
and each term is estimated in T~ a Stz\> by using 

\ u i u \\\m ^ II^iIIh/l-'II^iII// 1 ^ ** -1 ~ a ll«ll2t> 
l^ill^ < KlU*Kll^ <t- 2a ||u||| t , 

I«?Vt4||^ < hilU-l|Vn 2 || L2 < 5t' l ' a \\u\\ Zt , (5.15) 
\u\Vul\\ L 2 < ||^|| L2 ||Vu°|| L oo < St-^Mz,, 
\u\Vu\\\ LP < ||«i|U,||Vi4|| L2 < t- 2a ||n||| t , 

where 5 is the small factor coming from \\(p\\gi , and we choose p € (1,2) and g G (2, oo) 
such that 

a - k > 1 - l/q = 3/2 - 1/p. (5.16) 

Thus we obtain 

IIW(«)lls*4 < ^ a [6 + T- k \\u\\ 2t ] \\u\\ 2t , 

II WW - WMIU4 S + t ~ k (IMUt + IMM ll« - Hz-- 
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5.5. Iteration argument. We define an iteration sequence (zf k \ , U( k \ , Vr k \ ) for k 
and t > T > 1, by 

2(0) = v (o) =z° = e~ iHt <p, u (0 ) = u° = VV\ 
Z(k+1) = z° + Tri(u {k) ) + Dif(u (k) ) + ^sy(u°), 



!^(fc)l 2 



2 

We introduce the foilowing norm for z: 

\\ Z \\Z 2 := N f» ill + 2 L(3 rjl/2- 

Since /? < a, we have < T^H^Hg^. Using the Sobolev embedding H x ( 

(3 < a, we have for any k, j = 0, 1, 2 . . . , 

IKfc+i) -"cm-ijIU^, ~ \\v(k+i) - v u+i)\\ L % tT Hi + IK*+i) - n a+i)ll^ T L4 

~ \\ Z (k) - z (j)\\zl + lll^(fc)| 2 - Ki)| 2 |U« T £|5 

and then the quadratic part is estimated by using 1)5.9}) and a + k < 2/3, 

IH^Cfc)! 2 - \ u (j)\ 2 \\L« T Li ^ [||«(fc)||2r + IKj)lkr] ||«(fc) -U(j)\\z' T - 

For the first iteration, we have um — it(o) = — i- > |u°| 2 /2 S R, and so 

-«(0)II^ itL 4 < ||«(1) -«(0)IIl- t ^i £ III^TlU^ < T- K \\u°\\ 2 Zt . 

As for we use (|5.13p and l|5.17j) . deriving 

||^fc + i) - ^ + i)|| 5t 4 < T' a [5 + T~ K (1 + \\u {k) \\ ZT + \\u U) \\ Zt ) 2 ] 

x ||«(fc) -U(j)\\ z > T - 
For the first iteration, we apply (|5.13l) and (|4.2|1 . Then we get 

Hi)-mH £ \\Tri{u°)\\ T - aStz i T + \\Asy{u°)\\ Z 2 
<t- k [\\uYzt + , 

where we denote 

Hk:=Hifli+ E ii(e)~ 1/2 iei |fe| ^(e)iu 2 nLoc. 

0<|fc|<2 

Gathering the above estimates, we deduce that 

D k ■= \\z( k ) - Z(k-i)\\z2 + IKfc+l) ~ u {k)\\z' T , E k '■= \\ u (k)\\z T , 

Dk+i <(S + + E k+1 +E k + Sfe-i) 2 )^, (fc = 1,2,3... ), 

Di<T^[\\uYz T ^ + \\A\z T ) + \W\\lf], 
fe-i 

^ £ E^ + IKi)ll^> IK(i) IUt < IK IUt(i + ^" k II«°IU t ). 
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Hence for sufficiently small 5 and large T, (z^,u^) converges to some function (z,u) 
satisfying the equation (|1.15|) and 

\\ z -*(o) \\z* + + \\ z ~ z (i)\\t-<* 5*4 ^ D l- (5.27) 

The uniqueness for t > T is proved also by the above difference estimates. In addition, 
(|4~2^1 implies that 

t l - e {\\z'\\^ + \\z {x) - z° - z'\\ Ll )+f' 2 \\z'\\H e < M 2 M , (5.28) 

where we used the fact that zn\ — z° — z' = Asy'(u°) does not contain $o- By using this 
estimate together with (|5.22(l and the L p decay, we also obtain 

Mfrnfr Z \\u\\ 2 L * < t- 20 (D 1 + \\u°\\ Zt ) 2 , 

\H\h2s < IMlJa/Cl-J £ (\\z\\ L 2/(l-e) + \\U\\ 2 L 4) 2 (5.29) 

<t- £ (\\<p\\B 0/ + ||* + J D 2 + ||u°||| T ) 2 , 

for any small e > 0. In particular, z(T),v(T) G P s for < s < 1. The local uniqueness for 
(|1.17|) in this class is easily derived from the Strichartz, Sobolev and Holder inequalities. 



5.6. Global continuation. The final task is to extend our solution u to t < T. We can 
not apply the P 1 global wellposedness of u by 2 , since Asy(u°) barely falls out of L 2 . 
However the nonlinear energy is still finite, because u\, \ u\ 2 G L 2 , and the L 2 singular part 
at low frequency belongs to L°°. Hence we can apply the global existence results in [TO] , 
and we have only to see persistence of our function space, namely (z, v) G C(R; P e Pi P 1 ). 
The conserved energy can be written as 

E(a) . fEf^M^^fl^f^ ( , 30) 



where we denote p := u± + ^j-. jlUl Theorem 1.1] gives 

P eC(R;L 2 ), «£ C(R;(L°° + H^nH 1 ), 
u _ e iA(t-T) u(T) e c^H 1 ). 

Hence we have u G C(R; P s ) for < s < 1, which implies for < s < 1 that 

PC/ _1 |n| 2 G C(K;P 2s ). (5.32) 
On the other hand, by using the identities 



_> 



(5.31) 



\u\ 2 

Uz 1 = m + P^- = Pp + Qm, z 2 = u 2 , ( 5 - 33 ) 



we get 



z G C(K; (L°° + P 1 ) n P 1 ), Vp G C(K; L 2 + L 8 / 5 ), 



(5.34) 
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and the integral equation for z can be written as 

z =e- lH ^- T h(T) 

P rr, \ (5.35) 

- J e- lH{t - s) [2i{pui) + 2PC/ -1 V ■ (u 1 Vu 2 ) - 2U(pu 2 - u lU2 )]ds, 

where the nonhnearity is in C(R; H 1 + ff 1,4//3 ). Hence the Strichartz estimate implies that 

z{t) - e- iH(t - T) z(T) G C(M; F 1 ), (5.36) 

and therefore 

zeC(R;H s ), (5.37) 
for < s < 1. Combined with (|5.32j) . this implies that 

v = z - PU~ 1 \u\ 2 /2 £ C{R;H S ), (5.38) 

for < s < 1. □ 
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